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Dynamical decoupling is a coherent control technique where the intrinsic and extrinsic couplings 
of a quantum system are effectively averaged out by application of specially designed driving fields 
(refocusing pulse sequences). This entails pumping energy into the system, which can be especially 
dangerous when it has sharp spectral features like a cavity mode close to resonance. In this work 
we show that such an effect can be avoided with properly constructed refocusing sequences. To this 
end we construct the average Hamiltonian expansion for the system evolution operator associated 
with a single "soft" 7r-pulse. To second order in the pulse duration, we characterize a symmetric 
pulse shape by three parameters, two of which can be turned to zero by shaping. We express 
the effective Hamiltonians for several pulse sequences in terms of these parameters, and use the 
results to analyze the structure of error operators for controlled Jaynes-Cummings Hamiltonian. 
When errors are cancelled to second order, numerical simulations show excellent qubit fidelity with 
strongly-suppressed oscillator heating. 

PACS numbers: 03.67.Pp, 03.67.Lx, 82.56.Jn 



I. INTRODUCTION 

Quantum coherent control has found way into 
many applications, including nuclear magnetic resonance 
(NMR), quantum information processing (QIP), spin- 
tronics, atomic physics, etc. The simplest control tech- 
nique is dynamical decoupling (DD), also known as re- 
focusing. The goal of preserving coherence by averaging 
out the unwanted couplings is achieved most readily by 
runningprecisely designed sequences of uniformly-shaped 
pulses @, |,| . 

In a closed system, the corresponding performance 
can be analyzed in terms of the average Hamiltonian 
theory j|, | . To leading order, the evolution over the refo- 
cusing period r is indeed described by the time-averaged 
Hamiltonian of the system in the "rotating frame" de- 
fined by the control fields. Generally, the average Hamil- 
tonian is constructed as a series in powers of r. The 
number of the leading terms of this expansion that are 
exactly zero determines the order K of the refocusing se- 
quence. Larger K imply asymptotically more accurate 
refocusing, with error terms scaling to zero faster with 
decreasing r. 

For an open system, the dynamics associated with the 
bath degrees of freedom can be also averaged out, as long 
as they are sufficiently slow. With leading-order (K = 1) 
refocusing, the state decay processes are dramatically 
suppressed [| |, with a moderate decrease of the dephas- 
ing rate||, while with second-order refocusing (K = 2) 
both decay and dephasing are strongly suppressed ||. 

The decoherence analysis in Ref. || was based on 
the assumption of the low-frequency oscillator bath be- 
ing near thermal equilibrium. This assumption becomes 
questionable if the bath has sharp spectral features — e.g., 
if the controlled qubit system is coupled to a local high-Q 
oscillator. On the other hand, such a situation where the 
controlled system is coupled to an oscillator mode is quite 



common. This situation is realized in atomic physics, 
where the oscillator in question is the cavity mode, while 
the continuous- wave (CW) excitation is used to suppress 
the coupling^]]. In several quantum computer designs, 
nearly- linear oscillator modes are inherently present (e.g., 



mutual displacement in ion traps) 10 



QCs based on electrons on helium |15|, |lq , |17| , |lq|). 
nally, there are suggestions to include local high-Q oscil- 
lators in the QC designs to serve as "quantum memory" 
|l9| , "quantum information bus" [^0 , [21 , 2^] , or as a part 
of the measuring/control circuitry 



In this work we consider dynamical decoupling in a sys- 
tem where the spectral function of the oscillator bath has 
a sharp resonance. We include the resonant mode and 
the corresponding couplings in the system Hamiltonian, 
and consider the dynamics of the closed system driven 
by the refocusing pulses applied to the qubits only. We 
construct the average Hamiltonian for a situation where 
one of the qubits is driven by a single symmetrical one- 
dimensional 7r-pulse. To second order, the expansion is 
characterized by three parameters, two of which can be 
turned to zero by pulse shaping. An analysis of any refo- 
cusing sequence is then reduced to computing an ordered 
product of evolution operators for individual pulses. We 
illustrate the technique by analyzing the controlled dy- 
namics of a single qubit coupled to an oscillator. One of 
the analyzed sequences provides an order K = 2 qubit 
refocusing for any form of qubit-oscillator coupling. The 
simulations done for the Jaynes-Cummings Hamiltonian 
show excellent qubit fidelity with strongly-suppressed os- 
cillator heating, as long as the oscillator frequency bias 
exceeds the small coupling between the qubit and the os- 
cillator remaining in the effective Hamiltonian. We argue 
that results of Ref. || for corresponding open system re- 
main applicable as long as this renormalized coupling is 
small compared to the resonance width. 
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II. BACKGROUND 

A. Dynamical decoupling and effective 
Hamiltonian theory 

The main idea of dynamical decoupling is to drive the 
system in such a way as to average out the effect of 
unwanted Hamiltonian couplings. Obviously, this only 
works if the control fields are large compared with the 
other terms of the system Hamiltonian Hg. 

The easiest situation to analyze is where "hard" 6- 
f unction pulses are used. In this case the system Hamil- 
tonian can be ignored altogether during the action of the 
pulse. For a single qubit, a 7r-pulse along the cc-axis cor- 
responds to the evolution operator X = exp(— iirs x ) = 
—iux, where s x = Oxl^ is the spin-1/2 operator. The evo- 
lution operator for a sequence of such pulses interrupting 
periods of free evolution can be written as a product of 
the corresponding unitaries. For example, the standard 
spin echo|24[ sequence of a 7r-pulse and a negative 7r-pulse 
in the x direction followed by intervals of free evolution 
of equal duration r corresponds to the operator 



U 



X-T-X-T- 



(1) 



Such expressions are easily simplified using the corre- 
sponding matrix algebra. For the case of NMR, the sys- 
tem Hamiltonian is that of the chemical shift, 



H s = -Aa z , 



(2) 



it anticommutes with the pulse unitary X, thus 
e~ lHaT X = Xe +lHsT , and the two-pulse sequence (||) 
simplifies to the identity operator, 



U 



X-T-X-T- 



X e +iHsr e - iHsT X = XX 



1. 



(3) 



The simplicity of this formalism led to a number of strong 
mathematical results applicable to refocusing with ideal 
5-pulses. In particular, a succession of "concatenated" re- 
focusing sequences provide an excellent refocusing accu- 
racy which grows very rapidly with the number of pulses 
in a sequence p5|, p6| . 

In practice, however, the hard-pulse condition may be 
difficult to satisfy, and one has to account for the cor- 
rections associated with the action of the system Hamil- 
tonian Hg during the pulse. If we denote the control 
Hamiltonian Hc{t), the total Hamiltonian is 



H(t) = H c (t) + H s . 



(4) 



The simplest lst-order decomposition (e.g., see Ref. ||]) 
amounts to adjusting the intervals of free evolution before 
and after the pulse, 



-i{Hs + H c )r 



e — iHs Tl g - iHc T g - iHs T2 



(5) 



where Hq is assumed time- independent, with the precise 
value of Ti and r 2 computed in order to optimize the ac- 
curacy according to some fidelity measure. Superficially, 



any combination such that n + t 2 — t appears to provide 
equal accuracy to first order in r. In fact, the accuracy of 
the expansion relies on both H$t and Hqt being small; 
the results change non-trivially for finite-angle rotations. 

A more systematic way to analyse the effect of pulse 
shape is in terms of the average Hamiltonian theory [Q ||. 
This is equivalent to constructing the cumulant expan- 
sion of the evolution operator in powers of the system 
Hamiltonian in the interaction representation with re- 
spect to the control Hamiltonian Hq (t) which is treated 
exactly. 

For a system of qubits, the single-qubit control can be 
written most generally as 



Hc(t) 



(6) 



where er^, /i = x,y,z, are the usual Pauli matrices for 
the n-th qubit (spin). We assume that the qubit levels 
are nearly degenerate, or that Eq. ([|) is written in the 
rotating wave approximation with respect to the qubit 
working frequency, so that the system Hamiltonian docs 
not contain large terms. The control Hamiltonian (||) is a 
sum of single-qubit terms, and the zeroth order evolution 
operator which obeys the equation 



U {t) = -iH c (t)U (t), U o {0) - a, 



(7) 



can be constructed without much difficulty as a product 
of corresponding single-qubit operators. Then, the stan- 
dard prescription is to separate the fast dynamics due 
to control fields out of the evolution operator by the de- 
composition, U(t) = Uo(t)R(t), and write the equation 
of slow evolution for R(t), 

R(t) = -iH s {t)R(t), H s (t) = U^t) H s U (t). (8) 

The system Hamiltonian in the interaction representa- 
tion, Hs(t), is small at the scale of the refocusing period 
r, which makes the time-dependent perturbation theory 
(TDPT) expansion applicable. The Magnus expansion, 
the exponentiated version of TDPT, has somewhat better 
convergence properties. It is written in terms of cumu- 
lants Cfc, 



R(t) = exp(C 1 (t) + C 2 (t) + ■■■), 



Ci(t) 
C 2 (t) 



'o 

ft i-U 



I 



dt 2 / dti[H s {tx),H s {t 2 )l 



(9) 
(10) 

(11) 



Generally, the fc-th cumulant Ck contains a fc-fold inte- 
gration of the commutators of the rotating-frame Hamil- 
tonian Hs{U) at different time moments ti and has an 
order (tHs) k - 

Let us consider periodic dynamical decoupling with the 
period r, V£(t+T) = VJf(t). In addition, we request that 
the zeroth order evolution operator (Q) should also be 
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periodic, Uq(t) = 1 (zeroth-order refocusing condition). 
Then, the Hamiltonian Hs{t) is periodic, and the evolu- 
tion operator over a time interval commensurate with t, 
t„ = jit, can be factorized, 

R{tit) = [R{ T )] n = exp(-inrH ave ), (12) 

where the "average" Hamiltonian H avc is defined in terms 
of cumulants, 

-irH avc = d(T) + C 2 (t) + ■ ■ ■ . (13) 

The main advantage of the average Hamiltonian theory 
is the improved convergence at large time: regardless of 
the value of t, the TDPT series needs to be convergent 
only at \t\ < r. 



B. Order of a refocusing sequence 

For order- if refocusing, K > 1, we require additionally 
that the first K terms in the expansion of the effective 
Hamiltonian vanish, C\(t) = ... = Ck(t) — 0. For 
a closed system, a refocusing sequence of higher order 
generally offers better scaling of accuracy with the se- 
quence period. If we denote the maximum coupling in 
the system Hamiltonian as J (more precisely, J = |[i?s|[j 
a norm of the system Hamiltonian), the fc-th cumulant 
scales as Cf. ~ (Jr) k , while the corresponding term in the 
average Hamiltonian i/lvl ~ J fc r fc ~ 1 . Consequently, for 
order- if dynamical decoupling, the refocusing error de- 
fined as the norm of the deviation of the unitary operator, 
S rj = \\U(t) - 1||, scales as 5 V ~ (t/T)(Jr) K+1 , while the 
corresponding fidelity scales as 1 — F ~ Sfj oc t 2K . Thus, 
in the same system, a refocusing sequence of higher order 
can be run at a slower repetition rate. 

Additional advantage of the effective Hamiltonian the- 
ory comes from its locality. The first-order effective 
Hamiltonian (cumulant C\) contains only the qubit in- 
teractions present in the original Hamiltonian iis. The 
terms C2 contains connected pairs of such coupling terms. 
Generally, if we represent system Hamiltonian iig as a 
graph, with qubits as vertices and two-qubit couplings 
as corresponding edges, the expansion of the A:-th order 
cumulant Ct can be represented as connected subgraphs 
with up to k edges. As a result, lst-order refocusing con- 
dition, Ci(t) = 0, can be verified by analyzing clusters 
with up to two qubits; second-order refocusing, with the 
additional condition C*2(t) = 0, requires analysis of all 
clusters with up to three qubits, etc. 

This property, otherwise known as the cluster 
theorem p7|, allows one to design scalable refocusing se- 
quences whose refocusing order is independent of the sys- 
tem size. For a linear chain with nearest-neighbor in- 
teraction, one can achieve this by intermittently pulsing 
odd and even-numbered qubits. Several such particular 
sequences of order K = 2 and higher were demonstrated 
in Ref. M 



C. Sequence order and open-system refocusing 

Dynamical decoupling can be also effective against de- 
coherence due to low- frequency environmental modes [^9| . 
This can be understood by noticing that the driven evo- 
lution with period r = 27r/f2 shifts some of the system's 
spectral weight by the Floquet harmonics, to — > u) + nil. 
With the first-order average Hamiltonian for the closed 
system vanishing (if = 1 refocusing), the original spec- 
tral weight at n = disappears altogether, and the direct 
transitions with the bath degrees of freedom are also sup- 
pressed as long as f2 exceeds the bath cut-off frequency, 
^ ^ w c H 0- This corresponds to effective suppression 
of dissipative (Ti) processes. 

The full analysis of decoherence, including both dis- 
sipative and reactive (T2) processes, in the presence of 
order- if dynamical decoupling, if < 2, was done by one 
of the authors using the non-Markovian master equation 
in the rotating frame defined by the refocusing fields 
This involved a resummation of the series for the Laplace- 
transformed resolvent of the master equation near each 
Floquet harmonic, with subsequent summation of all har- 
monics. 

The results of Ref. M can be summarized as follows. 
With if > 1 refocusing, there are no direct transitions, 
which allows an additional expansion in powers of the 
small adiabaticity parameter, lu c /Q. In this situation the 
decoherence is dominated by reactive processes (dephas- 
ing, or phase diffusion). With K = 1, the bath corre- 
lators are modulated at frequency tt. This reduces the 
effective bath correlation time, and the phase diffusion 
rate is suppressed by a factor oc u> c /Cl. With if = 2 re- 
focusing, all 2nd-order terms involving correlators of the 
bath coupling at zero frequency, lu — 0, are cancelled. 
Generically, this leads to a suppression of the dephasing 
rate by an additional factor cx {uj c /Vt) 2 , while in some 
cases (including single-qubit refocusing) all terms of the 
expansion in powers of the small adiabaticity parameter 
(w c /f2) disappear. This causes an exponential suppres- 
sion of the dephasing rate, so that an excellent refocusing 
accuracy can be achieved with relatively slow refocusing, 

il > LU C . 



III. DYNAMICAL DECOUPLING OF A 
GENERIC ALLY-COUPLED QUBIT 

The analysis in Ref. |^ was done for a generic thermal 
bath with a featureless quasi-continuous spectrum char- 
acterized by the upper cut-off frequency uj c . The absence 
of sharp features justified an approximation where bath 
memory effects were essentially ignored at the scale of the 
decoherence time, although bath correlations at shorter 
times are crucial for describing the effects of dynamical 
decoupling. A sharp spectral feature, like a high-Q cav- 
ity mode, makes the starting point of the analysis^, the 
non-Markovian master equation involving only qubits, 
questionable, unless the decoherence time in the absence 
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of refocusing is long compared with the equilibration time 
of the high-Q mode. 

A. Model 

In this work we include any sharp quantum mode(s) 
into the "system" part of the Hamiltonian Hg, and con- 
sider the driven quantum dynamics of the resulting closed 
system. Compared with a system of qubits, a quantum 
oscillator admits a wider variety of linear or non-linear 
couplings. By this reason we begin with a model of a 
qubit with most general couplings, 

H s = a x A x + a y A y + a z A z + A , (14) 

where tr M are the qubit Pauli matrices and A u , v = 
0, x, y, z are the operators describing the degrees of free- 
dom of the rest of the system which commute with a jll 
[<7 M , A v \ = but not necessarily with each other. 

B. Pulse structure to linear order 

First, consider the qubit evolution driven by a one- 
dimensional pulse, 

H c = \a x V x {t), 0<t< Tp , (15) 

where the field V x (t) defines the pulse shape. The unitary 
evolution operator to zeroth order in Hg is simply 

U (t) = e— »*«/ 2 , cf>(t) = f dt 1 V x {t'). (16) 

Jo 

When acting on the spin operators, this is just a rotation, 
e.g., U (t)<jyU^(t) = (T y cos(p(t) +a z sin0(t). 

For inversion pulses with the net rotation angle 4>(t p ) = 
7r, with a symmetric shape, V x (r v —t) = V x (t), the aver- 
age of the cosine over pulse duration is zero by symmetry, 

1 f Tp 

(cos (j)(t)) p = — dt cos <f>(t) = 0. (17) 
t p Jo 

In such a case, the first two terms of the expansion, X = 
+ TpX^ + t^xW + . . ., of the unitary evolution 
operator X = U(t p ) in powers of pulse duration r p read: 

X^ = -ia x , (18) 
X (1) = -A x - a x A + is (a y A y + a z A z ) . (19) 

Here the dimensionless parameter 

s = (sin <p(t)) p (20) 

is the only one that characterizes the pulse shape in this 
order. 

We note that for 7r-pulses the sine of evolution angle 
is non-zero only over the duration of the pulse; the time 
intervals before the beginning and after the end of the 



pulse where V x (t) — do not contribute to the value 
of s. Because of that, this parameter can be viewed as 
a measure of the effective duration of the pulse. For 
example, while s = for an infinitely short <5-pulse, for a 
Gaussian pulse fl32| of width r, 

7T 1 / 2 2, 2 

G T {t + T p /2) = e-« /r , r«r p , (21) 

T 

this parameter is s ~ 1.5t/t p (see Tab. |). 

To create a pulse shape with effectively zero width to 
linear order, one can compensate for positive values of 
sin <j)(t) near the middle of the interval by making V x (t) 
somewhat negative near the beginning and the end of the 
interval. Such a first-order self-refocusing pulse was first 
suggested by Warren J|(| as a "Hermitian" shape, 

H T (t + r„/2) = G T (t + t„/2) 1 il^f ■ (22) 

where the precisely computed value 7 = 0.9609317217 is 
somewhat different from that in Ref. The correspond- 
ing fixed-length lst-order self-refocusing pulse shapes 
Sl, with all derivatives up to and including (2L — 1) st 
vanishing at the ends of the interval, L — 1,2, were 
constructed |2§|| in terms of their Fourier coefficients H|, 

V(t + t p /2) = A Q + ^2A m cos(mn p t), (23) 

m 

where the angular frequency Q p = 2~k/t p is related to 
the full pulse duration r p ; the coefficients A m are listed 
in Ref. Compared to Hermitian pulse shape, the main 
advantage of pulses Sl is their smaller power [the maxi- 
mum amplitude of the field V(t)]. 

For lst-order self-refocusing pulse shapes such as H T 
or Sl, the terms proportional to s in Eq. ( |l9|) disappear, 
and the expansion of the unitary operator simplifies to 

X s=0 = -ia x - t p (A x + <t x A ) + 0(t 2 p ). (24) 

C. Pulse structure to quadratic order 

The structure of the pulse to quadratic order is eas- 
ily computed as the next order in the TDPT. We have 
[cf. Eqs. ©, ©], 

X {2) = + Z -({A ,A x } + a x (A 2 + At)) 

+C([A), cr y A z - a z A y ] + i{A x , OyAy + a z A z }) 
s 

+ -({A Q , OyAy + a z A z } - i[A x ,a y A z - a z A y ]) 
+a(Al + A 2 z + ia x [A y ,A z ]), 

+ ^({A z ,A y }+ i <j x (A y ]+Al)), (25) 

where we parametrized the pulse shape in terms of two 
additional parameters 

a = (6{t-t')sm[(j)(t) -<f>(t')]) p , (26) 
C = (O(t-t') cos 0(t')) p , (27) 
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with the two-time averages over pulse duration 

r T p dt f Tp dt 1 
</(t,0)p= / - / -/(*,*')• (28) 
Jo "p jo T p 

We also use the notations [A, B] = AB — BA for the com- 
mutator and {A, B} — AB + AB for the anticommutator 
of two operators. 

The effect of the parameter a was studied previously 
in Ref. ^8|. The condition a = is necessary to obtain 
NMR-style one-dimensional second-order self-refocusing 
pulses. Indeed, the corresponding system Hamiltonian is 
just that of the chemical shift, Eg. (||), thus A x — A y — 
A = and A z = A/2 in Eq. Q. Then, the evolution 
operator to quadratic order in f = r p A/2 is simply 

X A = -ia x + is?a z + f 2 (a + ia x s 2 /2) + C(f 3 ); (29) 

the linear and quadratic corrections disappear entirely for 
second-order self-refocusing pulses such that a = s = 0. 
These are exactly the conditions used to design pulse 
shapes Ql (Ref. |2^). The index L = 1,2 denotes the 
parameter for the additional condition that the first 2L 
derivatives vanish at the ends of the interval, (0) = 
VM(t p ) =0,l = 0,...2L-l. 

The actual values of the parameters s, a, and C for 
several pulse shapes are listed in Tab. |. We note that for 
an ideal <5-pulse, the parameters s = a = 0, whereas £ = 
1/4 is not particularly small. For all "soft" pulse shapes 
listed, the values of £ are quite close to this value. The 
values of the parameter a are numerically small for all 
lst-order self-refocusing pulses with s = 0. The second- 
order self-refocusing pulses Q\, Q2 with s = a = may 
work as a perfect replacement of 5-pulses to second order 
accuracy pq. 



pulse 


s 


a/2 


C 


TTS(t - Tp/2) 








1/4 


Go. 05 


0.0744895 


0.0349708 


0.249476 


Go. 10 


0.148979 


0.0653938 


0.247905 


#0.05 





0.00153849 


0.249647 


#0.10 





0.00615393 


0.248589 







0.0332661 


0.238227 


S2 





0.0250328 


0.241377 


Qi 








0.239889 


Q2 








0.242205 



TABLE I: Parameters of several symmetric pulse shapes. The 
first line represents the "hard" ^-function pulse, Go. 05 de- 
notes the Gaussian |i2"| pulse with the width r = 0.05r p , see 
Eq. (H), #0.05 is the corresponding Hermitian[^o| pulse, see 
Eq. (p2), while Sl and Ql denote the 1st and 2nd-order self- 
refocusing pulses from Ref. |2g| ], with up to 2L — 1 st derivative 
vanishing at the ends of the interval of duration t p . 



IV. COMMON PULSE SEQUENCES. 

Transforming Eqs. (|l9|), ( p5| ) appropriately, we can 
now easily compute the result of application of any pulse 
sequence. In particular, the 7r-pulse X applied along the 
—x direction can be obtained from (—X) with the sub- 
stitution a —a. As a result, e.g., the expansion of the 
evolution opeator for the one-dimensional sequence XX 
[cf. Eq. (Q)] can be written as 

XX = 1 - 2iT p (A + a x A x - s(a y A z ~ a z A y )) 

-2t*(A + a x A x - s(<j v A z - a z A y )f + 0(r p 3 ), (30) 

or it can be re-exponentiated as evolution over time in- 
terval 2r p with the average Hamiltonian 

H^x = A + a x A x - s(a y A z - u z A y ) + 0(t 2 ). (31) 

We can attempt to correct for the terms proportional to 
s by using a longer sequence, e.g., XXXX. However, 
while the term is corrected in the leading-order effective 
Hamiltonian, we acquire a correction in the next order, 

ft xxxx = ®xA x — sTp {A x , <j y A y + <r z A z } 

+lST p [Ao, <J V A Z - a z A y ] + 0{t 2 ). (32) 

While the external Hamiltonian Aq cannot be averaged 
out by acting on the qubit, the term proportional to a x 
can be also suppressed with the help of two-dimensional 
sequences. The expression for the unitary Uy(t p ) = Y re- 
sulting from application of the pulse along the j/-direction 
can be easily obtained from Eqs. (|il|), ( p5[ ) by cyclic per- 
mutation of indices. Then, for example, the refocusing 
sequence 4p = 4p(xy) = XYXY corresponds to the 
effective Hamiltonian 

H 4p = A + - (<r x A z ~ a z A y ) ^-[A ,a x A x - a y A y ] 

-T P ^y(A 2 x + A 2 Z ) + Tp l -^[A z ,A y ] 

-T p ^^a z {A x ,A y } + 0(T 2 , S T p ), (33) 

where we dropped linear in r p terms proportional to s. 

The symmetric 8-pulse sequence 8s = YXYXXYXY 
produces the effective Hamiltonian 

(i 1 
-[A Z ,A X + A y ] + —(cr x A y — (T y A x ) 

+~a y {A x ,A y } + -a z {A y ,A z } 

+i-[A ,(T y A z +cr z A x + ^<J z A y - ^a x A z ]^j 

-^r(°y{A 2 x + A\) + i[A y ,A z }) + 0(r 2 ), (34) 

while the antisymmetric sequence 8a = YXYXXYXY 
corresponds to 

H 8a = A + ^{a x A z -a z A y ) + O(T 2 ). (35) 
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We note that in the latter case there is a leading-order 
term proportional to s but no terms in order r p ; this 
sequence produces 2nd-order refocusing already with 1st- 
order pulses. 



V. QUBIT IN A CAVITY 

To illustrate these results, consider a qubit placed in 
a lossless cavity with a single mode nearly-resonant with 
the qubit. We consider the simplest case where the sys- 
tem can be described by the Jaynes-Cummings Hamilto- 
nian, 



a z - g(tfa- + a+b), 



(36) 



where a± = (a x ± icr y ), and u) r and loq are the freque ncy 
biases for the cavity and the qubit respectively. Eq. ( J3q ) 
can be also written in the form (|l4| ) with A x — —g(b + 
6 t )/2, A y = ig{tf - b)/2, A z = uj , and A = uj r tfb. We 
concentrate on the special case lu — which corresponds 
to working in the "rotating frame," with the control fields 
[Eq. (||)] applied on resonance with the qubit. Note that 
we assume the control fields to be applied directly at the 
qubit and not at the oscillator (23); for a linear oscillator 
the corresponding compensation can be achieved by a 
spectral filter. 



A. Sequence 4p 

When the sequence 4p is used with a Gaussian or other 
pulse shape with s ^ 0, the effective Hamiltonian to lead- 
ing order can be written as 



isg 



b) + 0{r p ). (37) 



The original exchange-like oscillator coupling in Eq. 
is replaced by the qubit phase coupling in Eq. ([57]). While 
the coupling magnitude is reduced by a small factor oc s, 
it does not go down with more frequent pulse application. 
The same holds true if the oscillator is replaced by an 
auxiliary qubit (left panel in Fig. [I]), or if the pulses are 
applied in the x-z direction [dashed lines in Fig. [l]; the 
effective Hamiltonian is given by Eq. (^) below]. One 
can see from Fig. [j] that the refocusing accuracy is more 
or less similar for all these cases. 

When the sequence 4p is applied with lst-order self- 
refocusing pulses (Hermitian or S*l), the first-order terms 
in the effective Hamiltonian are gone. This leaves terms 
linear in t p as the leading-order correction 



H, 



4p 



Wr6 t 6 _<Z£|*: . x ( 6 t_ 6 ) 



. 1 + 4C 



r P g 2 a z (b 2 - h.c.) + 0(s, ar p , r p 2 ) (38) 



where we set s = and dropped the relatively small 
term oc ar p . We note that the first two terms in Eq. (|38|) 




FIG. 1: (color online) Evolution of the qubit fidelity and the 
number of quanta at the oscillator (n) (minimized and max- 
imized over initial qubit states, respectively) at the end of 
the refocusing interval. Length-4 sequence 4p with Gaussian 
pulses Goio applied on resonance with the qubit, Wo = in 
Eq. (^6|). Left panels: oscillator restricted to n = 0, 1 states, 
which makes it effectively a qubit. Right panels: oscillator 
restricted to n < 8 levels, which over the simulation time 
is equivalent to infinity. Bottom panels correspond to oscil- 
lator in resonance with the qubit, top panels correspond to 
oscillator frequency bias u> r = 0.117(27r/r p ). Red solid lines 
correspond to control pulses applied along X and Y axes [se- 
quence 4p(xy)]; blue dashed lines correspond to control pulses 
applied along X and Z axes [sequence 4p(xz)]. 
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FIG. 2: (color online) As in Fig. [I] but with lst-order self- 
refocusing pulses Si. The corresponding curves for pulses Qi 
(not shown) are virtually identical with linear scale. 



disappear with uj r = 0. The third term is suppressed 
when the oscillator is replaced by an auxiliary qubit 
(b — > T-, — > r + , &t& — > t z /2). The same hap- 
pens if the sequence is applied in the x-z direction [se- 
quence 4p(a;.z) = XZXZ], with the corresponding effec- 
tive Hamiltonian 



H, 



(xz) 
4p 



+i M± CTx ( 6 t _ 6) + o(sr p , ar p , r p 2 ). (39) 



The net effect is that the refocusing accuracy is dramati- 
cally improved when the oscillator is replaced by an aux- 
iliary qubit, or when the four-pulse sequence is applied 
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FIG. 3: (color online) As in Fig. ^ but for the sequence 8a. 
The term oc s in Eq. © produces large errors comparable 
for those for 4p sequence in Fig. [j]. This indicates that this 
sequence is less stable to pulse shape errors as, e.g., the sym- 
metric sequence 8s. 



in the x-z plane, with the additional mode being in reso- 
nance with the qubit. In such cases, the only remaining 
linear in t p terms in the effective Hamiltonian are those 
proportionaly to the small parameter a, and the refocus- 
ing gets much more accurate, see Fig. ^. 



B. Sequence 8a 

The antisymmetric eight-pulse sequence 8a produces 
the effective Hamiltonian 

H Sa = w r tfb + ^(&t - b) + 0(r 2 p ). (40) 

The refocusing is lst-order with Gaussian pulses (Fig. ||), 
with the errors comparable to those of the sequence 4p 
(cf. Fig. |l]). With self-refocusing pulses Si [Fig. Q or Qi 
(not shown) the sequence produces 2nd-order refocusing. 
With all 1st and 2nd-order error terms cancelled, the 
refocusing accuracy is improved substantially. 
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FIG. 5: (color online) As in Fig. [l] but for the sequence 8s. 

C. Sequence 8s 

The symmetric eight-pulse sequence 8s produces the 
effective Hamiltonian 



u r tfb-^^a y (tf + b) 2 + 0(sT p , rl) (41) 



which corresponds to Eq. ( pq ) with the larger terms al- 
ready suppressed. The refocusing is first-order with ci- 
ther Gaussian [Fig. || or lst-order pulses [see Fig. || with 
pulses Si], and second-order with second-order pulses Ql 
(not shown). With the linear scale of our plots, there is 
only a slight difference between zeroth and first order 
pulses (due to higher-order terms), and no visible differ- 
ence between pulses Si and Qi. 

We note that the performance of this sequence is very 
close to that of the antisymmetric sequence 8a. Never- 
theless, the symmetric sequence 8s provides better sta- 
bility with respect to the pulse shape errors, and it is also 
expected to result in better visibility (initial decoherence 
in Ref. B) when used in an open system. 



D. Discussion 
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FIG. 4: (color online) As in Fig. [j] but for the sequence 8a 
with pulses Si. The corresponding curves for pulses Qi (not 
shown) are virtually identical with linear scale. 



Our simulations show that with properly designed re- 
focusing, to an excellent accuracy, the quantum oscilla- 
tor coupled to the qubit remains in the ground state, 
while the qubit fidelity remains very close to unity. Nu- 
merical results agree with the predictions from the an- 
alytically computed effective Hamiltonians which show 
greately reduced coupling between the qubit and the os- 
cillator. Then, in the closed system, a very small fre- 
quency bias (of order of renormalized coupling value) be- 
comes sufficient to effectively disconnect the oscillator 
mode and protect the qubit coherence. Similar effect can 
be achieved when the finite width T of the resonance is 
taken into account: for refocusing to be effective, T must 
be large compared with the renormalized qubit coupling. 

In other words, refocusing is effective as long as the 
renormalized qubit coupling with the resonant mode is 
small on the scale of the resonance width. This gives the 
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FIG. 6: (color online) As in Fig. [l] but for the sequence 8s 
with pulses Si. Formally, the sequence is of the same order as 
with Gaussian pulses, see Fig. [|; the noticeable differences are 
due to error terms of higher order which we dropped in our 
calculations. The off-resonance refocusing is excellent (the 
fidelity and (N) curves in the top right panel run along the 
corresponding axes), while the on-resonance performance is 
also good. The plots with second-order pulses Qi (not shown) 
look almost identical. Note also that with the linear scale of 
these plots, the curves here are almost indistinguishable from 
those in Fig. [l|. 



modified, more broadely applicable, criterion for appli- 
cability of the results of Ref. In particular, we ex- 
pect the second-order sequence 8s with 2nd-order pulses 
Q„ to provide an excellent refocusing accuracy to system 
Hamiltonian of Jaynes-Cummings form ( |36| ) also in the 
presence of a thermal bath, as long as the refocusing rate 
is sufficiently high. In addition to the condition on renor- 
malized value of the high-Q oscillator mode coupling, the 
refocusing period r must be below the threshold value of 
order of the inverse bath cut-off frequency refocusing fre- 



quency, lu c , see Ref. 



VI. CONCLUSIONS. 

In this work we analyzed the performance of soft- 
pulse dynamical decoupling in the presence of a sharp 
resonance mode. Because of the associated memory ef- 
fects, such a problem cannot be addressed by considering 
master-equation dynamics for the qubit-system density 
matrix alone. Instead, we included the resonant mode 
and the corresponding couplings in the system Hamilto- 
nian, and considered the dynamics of the resulting closed 



quantum system driven by a sequence of soft or hard re- 
focusing pulses applied to the qubit. The analysis was 
done in terms of the effective Hamiltonian theory which 
describes the evolution of the system "stroboscobically" 
at the time moments commensurate with the refocusing 
period. 

In fact, to make our results applicable to a large num- 
ber of possible coupling terms between the qubit and the 
oscillator, we solved a more general problem of an arbi- 
trary coupled [see Eq. (|l4|)] controlled qubit. The main 
result of this work is the expansion of the unitary evo- 
lution operator (|l9|), ( p5| ) and the classification of the 
corresponding parameters in Tab. [|. This allows an ex- 
plicit computation of the error operators associated with 
refocusing in systems of arbitrary complexity. 

We also computed the effective Hamiltonians for sev- 
eral single-qubit refocusing sequences. To quadratic or- 
der, all coupling terms are cancelled if the length-8 se- 
quences 8a [Eq. (|35|)] or 8s [Eq. are used with 1st 
or 2nd-order self-refocusing pulses respectively. We illus- 
trated the general analytical results on the specific ex- 
ample of the driven Jaynes-Cummings Hamiltonian. The 
results of simulations agree with the predictions based on 
the analytically computed second-order effective Hamil- 
tonians, although in some cases the effects of higher-order 
terms not included in the calculation are noticeable. 

For robust single-qubit refocusing, we recommend the 
symmetric 8-pulse sequence 8s applied with second-order 
self-refocusing pulses. This sequence provides excellent 
refocusing accuracy for any form of the coupling of the 
qubit with outside world, and it is stable with respect 
to pulse shape errors. As the second-order symmetric 
sequence, it should also work well for open systems, as 
long as the environment is slow on the scale of the refo- 
cusing rateQ. The thermal bath is expected to remain 
close to equilibrium, and refocusing to perform well, as 
long as sharp features in the spectral function are wide on 
the scale of the renormalized value of the corresponding 
couplings. 
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